Exact results for the criticality of quench dynamics in quantum Ising models 
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Based on the obtained exact results we systematically study the quench dynamics of a one- 
dimensional spin- 1/2 transverse field Ising model with zero- and finite-temperature initial states. We 
focus on the magnetization of the system after a sudden change of the external field and a coherent 
time-evolution process. With a zero-temperature initial state, the quench magnetic susceptibility 
as a function of the initial field strength exhibits strongly similar scaling behaviors to those of the 
static magnetic susceptibility, and the quench magnetic susceptibility as a function of the final field 
strength shows a discontinuity at the quantum critical point. This discontinuity remains robust 
and always occurs at the quantum critical point even for the case of finite-temperature initial 
systems, which indicates a great advantage of employing quench dynamics to study quantum phase 
transitions. 

PACS numbers: 64.60.A-, 03.75.Lm, 03.75.Kk, 42.50.Pq 



I. INTRODUCTION 

The study of quantum phase transitions (QPTs) is 
a fascinating topic in condensed matter physics and 
quantum information science. QPTs distinctively from 
temperature-driven critical phenomena occur due to the 
competition between different parameters describing the 
interactions of the system, and QPTs occur only at zero 
temperature In the second-order QPT, the ground 
state undergoes qualitative changes when an external 
parameter passes through the quantum critical points 
(QCPs). The QCPs are characterized by the diver- 
gence in the correlation length, which leads to the crit- 
ical scalingbehaviors governed by a class of universal 
exponents [2|]. In principle, experimental observations of 
QPTs could be achieved at zero temperature How- 
ever, in practice, it is difficult to realize since cooling 
matter to zero temperature is impossible in any experi- 
ment. Hence, recently, the finite-temperature properties 
of QPT systems begin to attract more attention Q. In 
their results, the quantum criticality can persist up to a 
surprisingly high temperature. However, the critical be- 
havior at finite temperature is not exact but a remanent 
of that at zero temperature. Some divergent physical 
quantities at the QCP, e.g., the magnetic susceptibility in 
the one-dimensional spin- 1/2 transverse field Ising model 
(TFIM), become convergent and shifted away from the 
QCP in the parameter space at finite temperature [2(. 

On the other hand, besides the ground state, the dy- 
namic properties of QPT systems also cause a lot of in- 
terests 0,1, i 0, 1 I Efl M M E3- Actually, the 
zero-temperature static properties are naturally and in- 
timately linked to the dynamic process, in which many 
excited states are involved. As we know, the scaling ex- 
ponents are dependent on the effective dimensionality, 
which is the sum of the dimension and the dynamic ex- 
ponent 0, Q. For the adiabatic approach of the QPT 



system to the QCP, the variation of parameter needs to 
slow down to infinitesimal 7]. So the ground state at 
the QCP could not be achieved adiabatically. Differently, 
quench dynamics can pass through the QCP without any 
restriction. In experiments, the quench dynamics of the 
Bose-Hubbard model have been performed with ultracold 
bosonic atoms in optical lattices @. Quench dynamics 
and their critical behaviors in the TFIM are also investi- 
gated theoretically [1(1 (HI- However, none of them dis- 
cussed the effect of temperature on these critical behav- 
iors. The quench dynamics relate to all of the eigenstates 
more than just the ground state, so the quench dynamic 
properties at the QCP should not be as sensitive to tem- 
perature as those of the ground state. 

With these motivations, in this paper, we study the 
quench dynamics of the TFIM with zero-temperature and 
finite-temperature initial states. Initially, the system is 
modulated to be with a transverse field strength Xi and 
temperature T. The field strength is suddenly changed to 
A / and the system begins to evolve coherently. Here, we 
use the magnetization per spin in the transverse field di- 
rection m(Aj, Xf,T;t) to characterize the state at time t. 
After a long enough time evolution, the off-diagonal con- 
tributions to m(Aj, A/,T; t) are cancelled with each other 
flij , and m(Aj ,Xf,T;t) achieves a steady value finally. 
The final asymptotic magnetization is 



j(X h X f ,T) = m(Aj, \ f , T; oo), 



(1) 
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which is called the quench magnetization in this paper. 
To make the above formula more easily tractable math- 
ematically, we write it as [3] 

r dt 

m q (X i ,X f ,T)= lim / -m(X h X f ,T;t). (2) 

We will show that m q as well as its derivatives, the 
magnetic susceptibilities Xii^Hi ^f>T) = dm q /dXi and 
X/(Ai, A/, T) = dm q /dXf, exhibit critical behaviors when 
Aj or Xf passes through the QCP. For a zero-temperature 
initial state, the susceptibility Xi diverges logarithmically 
and exhibits scaling behaviors when Xi is in the vicinity 
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of the QCP. For an initial state with temperature T, the 
susceptibility Xf experiences a jump Sxf = tanh(Aj/2T) 
when Xf is at the QCP, where A; is the energy gap be- 
tween the ground and first excited states of the initial 
system. The jump 8\f achieves maximum when T turns 
to zero. Since this jump is not sensitive to T and the 
operations are not restricted by adiabatic conditions, the 
quench dynamic process has the advantage to study crit- 
ical behaviors of QPT systems without the rigorous re- 
striction of zero temperature. 

II. QUENCH MAGNETIZATION 

The TFIM is a famous model for studying the second- 
order QPTs. It is exactly solvable and useful for verifying 
many new concepts and methods. The Hamiltonian of 
the TFIM is 

N 

ff(A) = -5>J^ 1 + AaJ), (3) 

where a a is the Pauli matrix (a = x, y, z) and N is 
the number of sites. The QCP of this system is at 
A = 1. Initially, the strength of the transverse field is 
Aj and the system is in a thermal state p{\,T;Q) — 
Z~ x exp[— H(Xi)/T] with temperature T, where Z = 
Trexp[— H(Xi)/T]. Then the transverse field is suddenly 
changed to Xf and the system begins to evolve coher- 
ently driven by H(Xf). The magnetization per spin in 
the transverse field direction is 

m{X i ,X f ,T;t) = TrX x p(X i ,T;t), (4) 
where Y* x = l/NJ2f =1 °j and p{X l ,T;t) = 

e -iH(\ f )t p ^. )T . j e iH(\ f )t is the state at time t 

After a long enough time evolution, m approaches the 
steady quench magnetization 

m q (Xi,Xf,T) = Trm 9 (A / ) j o(A i , T; 0), (5) 

where the quench magnetization operator fh q is 

m 9 (x f ) = E K A /)> M A /)I s * K A /)> M A /)I • ( 6 ) 

n 

Here, {|n(A))} is the complete set of eigenstates of H(X), 
i.e., H(X) |n(A)) = E n (X) \n(X)}, and n = corresponds 
to the ground state. 

Using the exact solution of TFIM [l], [f| [l5[ and taking 
the thermodynamic limit N — > oo, we get the quench 
magnetization as 

r 27! dk 

m q (X u X f ,T) = J — tanh[e(Ai, *)/2T] (7) 

x cos#(A/, k) cos[0(A/, k) — 9(Xi, k)], 

where the function 9(X,k) is defined as cos#(A, k) = 
2(A — cosfc)/e, sin0(A, k) = 2sinfc/e and e(A, k) = 
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FIG. 1: (Color online) Scaling behaviors of the susceptibility 
Xi(Xi, Xf ,0) with Xf = 0, 0.6, 1.5 and oo. (a) and (b) are 
plots of Xi(^i; Xf, 0) in the thermodynamic limit, (c) and (d) 
are plots of Xi(Xm, Xf ,0) and In |A m — 1| as functions of IniV, 
where A m is the pseudo-critical point. The slopes of the solid 
lines are — 7r -1 in (b), 7r _1 in (c), and —1.9 in (d). 

2V1 + A 2 - 2A cosfc. The quench magnetization is de- 
pendent on two systems, the initial system with the 
transverse field strength A^ and the final system with A/. 
In the following, we will consider two cases, A^ ~ 1 and 
Xf ~ 1, respectively. 

III. SCALING BEHAVIORS 

We firstly consider the case of zero-temperature initial 
state with A l - 1, where p(Ai,0,0) = |0(A^)) (0(A l )|. It 
should be noted that in the thermodynamic limit, the 
ground state of H(Xi) is two- fold degenerate in the re- 
gion Aj < 1. Fortunately, any combination of these two 
degenerate states as an initial state will give the same 
result. From Eq. (|5|) we have 

m q (X l , A/, 0) = (0(A l )|m,(A / )|0(A l )), (8) 

which writes m q (Xi, A/, 0) as an expectation value of the 
operator rh q (Xf) in the ground state |0(Ai)) and pre- 
sumably implies that a sudden change of the ground 
state |0(Ai)) around Xi — 1 will lead to some critical 
behaviors of m q (Xi, Xf,Q). Actually, in the strong-field 
limit Xf = oo, the transverse field term is dominant, 
which keeps the conservation of the expectation value 
of the total spin component in the transverse field di- 
rection during the quench process. Consequently, we 
have m q {Xi, oo, 0) = m q (Xi,Xi 1 0) and Xi(Ai,oo,0) = 
Xi(Xi,Xi,0), where m q {X l ,X i ,Q) and ^(A^A^O) are 
the static magnetization and susceptibility, respectively. 
Xi(Ai,oo,0) should share the same critical behaviors as 
those of Xi(Xi, Ai,0), the second order derivative of the 
ground-state energy density. So Xi(Ai,oo,0) should di- 
verge at Xi — 1 and obey scaling behaviors at A,- *~ 1. 
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Next we will show that a similar conclusion could be 
obtained for finite A/. Define A, = minjAjjA^ 1 } and 
Af = min{A/, AJ 1 }. When Aj > A/, from Eq. © we 
have [IB] 

/" Al dx ( v [k~\ 
m q (X i ,X f ,0) = ^ —l- + w ^—\ +w , (9) 

where v = sgn(A/ — x)c(Xf,x)[s(Xf,x)s(Xi,x) 
— c(Xf,x)c(Xi,x)], c(X,x) — [2A — (x + x~ 1 )]/d, s(X,x) 
= (x - x-^/d, d(X,x) = 2y/\l + A 2 - X{x + x^)\ 
and w — (Xi + A/)/ (27rA/\/ Aj Aj) . Accordingly, the 
susceptibility is 

^ ^ f Ai dx s(Xi,x) u , . . „. 

MX i ,Xf,0) = j Q -^ ) - + 0(X i ,X f ), (10) 

where u = sgn(A/ — x)c(Xf,x)[s(Xf,x)c(Xi,x) 
— c(Xf,x)s(Xi,x)] and C(Aj,A/) is convergent at 
Aj = 1. In the case of Aj ~ 1, i.e., Aj ~ 1, the above 
integral can be reduced as f Q *(...)ete ~ f.*_ e (...)dx with 
£ -C 1. Via straightforward calculations, an asymptotic 
behavior of Xi(Aj, A/, 0) around Aj = 1 is obtained as 

Xl (Aj, X f , 0) « -i In |Aj - 1| + K^Xf), (11) 

7T 

where K\(Xf) is a A/-dependent constant. In Fig. Q] (a) 
and (b), Xj(Aj, A/, 0) are plotted as functions of Aj and 
ln|Aj — 1| with Xf — 0, 0.6, 1.5 and oo. 

According to the finite size scaling ansatz pjj, the 
above critical behavior can be extracted from finite sam- 
ples, which is important for quantum simulations. Nu- 
merical simulations for finite systems show that the sus- 
ceptibility reaches the maximum Xj(A m , A/,0) at the 
pseudo-critical point A TO . In Fig. Q] (c), we plot 
Xi(A m ,A/,0) as a function of lnV with A/ = 0, 0.6, 1.5 
and oo. As expected, Xi{^m, A/,0) diverges logarithmi- 
cally as 

Xj(A m ,A J ,0)«-lniV + K 2 (A / ) ) (12) 

where K2(Xf) is another A/-dependent constant. As the 
lattice size approaches infinite, the pseudo-critical point 
A m tends to the QCP as |A m — 1| oc N^ 1 - 9 , which is shown 
in Fig. Q] (d) . According to the scaling ansatz in the case 
of logarithmic divergence [l7j , the ratio between the two 
prefactors of the logarithm in Eq. (|11[) and (fT2"|) is the 
exponent v that governs the divergence of the correlation 
length. As expressed in Eq. (|11[) and (fT2"|) . numerical 
calculations give v = 1, which accords with the result 
obtained from the exact solution of the TFIM (l~7j . 

When the initial system is at finite temperature, the 
correlation length of p(Aj, T, 0) is always convergent and 
all the scaling behaviors will vanish. We do not discuss 
this case in detail. 



IV. DISCONTINUITY OF THE QUENCH 
SUSCEPTIBILITY 

Now we focus on m g (Aj,A/,T) and %/(Aj, A/, T) at 
Xf ~ 1. Mathematically, under the condition A/ > Aj, 
Eq. (JT]) becomes 

|1- A 2 ,| 

m q (Xi,X f ,T) = J ' tanh8[c(Aj,A/) (13) 

-s(A l ,A / )] + A{Xi,X f ,T), 

where O = d(Xi, A/)/2T and A(Xi, Xj, T) is an analytical 
function when Aj is away from 1 [16j |. Obviously, the 
factor |1 — A* | in Eq. (fl3|) causes a sudden change of 
m q (Xi, Xf,T) at Xf = 1. For two extreme cases, Aj = 
and oo, Eq. (7} is completely integrable and could be 
written as 

m g (0, Xf,T) = tanh(l/T) { ^ J \ (14) 

and 

"» 4 (cx>,A /> r) = | 1 _ 1/(2A 2 )) A/>1 . (15) 

Accordingly, such a sudden change leads to a disconti- 
nuity of ^f,T) at the QCP. The expression is 

X/(Aj,A/,T) = sgnCAj-^^I^KAj.A/) (16) 

-s(Aj, A/)] + B(Aj, A/, T), 

where B(Xi,Xf,T) is continuous at A/ ~ 1. The sgn 
function leads to a jump of x/(^ij -Vj ^) with magnitude 
(5%/ = tanh(Aj/2T), where Aj = 2|Aj — 1| is the energy 
gap above the ground state for the initial Hamiltonian 
H(Xi). Remarkably, the jump 5\f always occurs at the 
QCP even for a finite-temperature initial state. When 
T <C Aj, Sxf decays slowly as 1 — 2exp(— Aj/T), and is 
almost a constant within the range T < O.lAj. In this 
sense, the lower-temperature samples share the same fea- 
ture as that of the zero-temperature sample, which is cru- 
cial for the experimental detection of critical behaviors of 
QPT systems, since cooling matter to zero temperature 
is impossible in any experiment. 

To exhibit the above critical behaviors, in Fig. [2l 
m g (Aj,A/,0) and x/(^i>^/)0) with Aj = 0, 0.5, 2 and 
oo are plotted. It is shown that m 9 (Aj, A/,0) has a sud- 
den change for systems with a wide range of Aj. For 
Xf < 1, m g (Aj,A/,0) increases as A/ increases, except 
that when Aj = oo, m g (Aj,A/,Q) is always equal to 0.5. 
In contrast, for A/ > 1, m g (Aj,A/,0) increases for Aj > 1 
and decreases for Aj < 1. Accordingly, x/(Aj, A/,0) has a 
discontinuity at A/ = 1 in the following manner: it drops 
down for Aj < 1 and jumps up for Aj > 1. 

It is well known that the ground state of the sys- 
tem H(Xf) experiences a sudden change when A/ passes 
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FIG. 2: (Color online) Plots of m q (Xi,Xf, 0) and x/O, A/, 0) 
with \i = 0, 0.5, 2 and oo. m q (Xi, A/,0) has a sudden change 
at Xf = 1, while X/(Ai, A/,0) has a discontinuity at A/ = 1. 
When A/ passes through 1, X/(Ai, A/,0) drops down for A; < 1 
and jumps up for A^ > 1. 
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FIG. 3: (Color online) The time evolution of m(Ai,A/,T, t) 
with Xi = 0, 2; T/Aj = 0.1, 1 and A/ = 0.8, 1, 1.2. 
m(Ai, Xf,T,t) oscillates initially and becomes steady at t = 15 
approximately. 



through the QCP. A natural question is whether the sud- 
den change of the ground state causes the critical behav- 
iors of quench quantities directly. Consider the simplest 
case with zero-temperature initial state, where the jump 
Sxf gets the maximum 1. From Eq. ([S]), the contribution 
of the final ground state |0(A/)) to the quench quantity is 
proportional to the fidelity of two ground states [HI, |Ti| 



KOCA.,) |0(A y )>| 2 = ncos 



(17) 



always vanishing for finite |A/ — A,;| in the thermodynamic 
limit. Thus the critical behaviors of quench quantities 
are not direct consequences of the sudden charge of the 
ground state, which indicates that the excited states also 
experience drastic changes at the QCP. 

V. DISCUSSION AND CONCLUSION 



In above analyses, we neglect the interaction between 
the sample of TFIM and its environment. Generally, 
such an interaction will induce decoherence of the system. 
Therefore, our results are valid only when the quench re- 
laxation time tq is short enough compared to the deco- 
herence time. In order to estimate the order of tq, nu- 
merical simulations for m(A i , Xf,T;t) are performed. In 
Fig. El we plot the results with X t = 0, 2, T/A, = 0.1, 1 
and Xf = 0.8, 1, 1.2 as examples. Remarkably, tq is not 
sensitive to T, Xi and A/, and tq is smaller than 15J -1 , 
where J is the Ising coupling strength and in this paper 
,/ = 1. In a spin network for quantum information pro- 
cessing (QIP), the shortest period of time for an opera- 
tion between two neighbor qubits via a natural dynamics 
is (2 J) -1 . Therefore the realization of a complete quench 
process is not a difficult task for a spin network which is 
utilized for QIP in practice. 

In summary we found that the quench magnetic sus- 
ceptibility as a function of the initial field strength ex- 
hibits strongly similar scaling behaviors to those of the 
adiabatic process, and the quench magnetic susceptibil- 
ity as a function of the final field strength shows a dis- 
continuity at the QCP, which remains robust even when 
the initial system is at finite temperature. This observa- 
tion is useful for understanding QPTs and studying the 
properties of the QPT systems. Moreover, it gives us a 
new approach to observe the QCP and critical behaviors 
using low-temperature samples experimentally, avoiding 
the rigorous restriction of zero temperature. Because of 
the universality principle, the critical behaviors only de- 
pend on the dimension and the breaking symmetry, so 
the obtained results are heuristic and may be extended 
to other QPT models. 



Straightforward calculations show that |(0(Ai) |0(A/))| is 
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